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Abstract 

Using the matrix product formalism we formulate a natural p-species generalization of 
the asymmetric simple exclusion process. In this model particles hop with their own 
specific rate and fast particles can overtake slow ones with a rate equal to their relative 
speed. We obtain the algebraic structure and study the properties of the representations 
in detail. The uncorrelated steady state for the open system is obtained and in the ( 
p — > oo) limit, the dependence of its characteristics on the distribution of velocities is 
determined. It is shown that when the total arrival rate of particles exceeds a certain 
value, the density of the slowest particles rises abroptly. 
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1 Introduction 



The Matrix Product Ansatz (MPA), first introduced in the work of Hakim and Nadal 
[1], become popular in the study of one dimensional stochastic exclusion models after 
Derrida et al [2] applied this technique to the simplest such model,namely the totally 
asymmetric simple exclusion process or ASEP[3]. For a recent review see [4]. They were 
able to find nontrivial representations for the relevant algebra and calculate many phys- 
ical properties of this process. Since then this technique has been applied to many other 
interesting stochastic systems, both with periodic and with open boundary conditions. 
Janowsky and Lebowitz [5] have considered the case where there is a single impurity (a 
fixed blockage) on the lattice, where hopping rates are reduced, and have found that this 
single blockage has global effects, due to long range correlations. The case of a moving 
blockage and the formation of shocks has also been addressed in several works [6-7]. A 
model consisting of two species moving in opposite directions have been studied in Evans 
et al [8] where spontenous symmetry breaking has been observed. Finally MPA has been 
formulated for models with quasi parallel updatings [9-11]. 

In most of these works , a simplified model of traffic flow [12] is cited as a physical appli- 
cation, specially when impurities are present which when encountered decrease the rate 
of hopping of other particles [13]. Evans [14] has applied MPA to the case of several kinds 
of particles hopping with their own rates on a ring and has found that above a critical 
density a macroscopic number of holes are condensed in front of the slowest particle. 
In the work of Evans, however, particles can not overtake, and the order of particles is 
unchanged during the process. 

Turning to the mathematical aspect of the problem, a matrix product state can be 
understood as a generalization of ordinary factorizable states with a product measure in 
which numbers are replaced with non-commuting objects. This allows the original prob- 
lem to be reformulated in terms of algebraic relations. The main advantage of this 
thecnique is that once a nontrivial representation is found, all the physical quantities 
like the currents , densities and correlation functions can in principle be obtained from 
the representation. This is of course not a simple task and one should overcome many 
difficulties of combinatorial type in calculating explicit form of the matrix elements of 
product of operators. 

Whether or not the zero energy state of a Hamiltonian with nearest neighbor interaction 
can always be formulated as a matrix product state has been answered in the affirmative 
by Sandow [15]. Thus to any Hamiltonian of the above type, corresponds an algebra 
in the MPA formalism. However finding nontrivial representions of this algebra may 
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be harder than the original problem itself. Therefore some authors [16,17] have tried 
to reverse the problem, thats to begin with a quadratic algebra,with hopefully simple 
representations and then to find the relevant stochastic process. In [17] the algebras 
considered are appearantly related to the partially assymetric case. This restriction has 
a two-fold disadvantage . First they obtain the concrete form of the algebraic relations 
only in cases for low number of species, and second they do not obtain [17] " a general 
framework for arbitrary N ( number of species ), which when specialized to N=2 gives 
the know results." 

Owr aim in this paper is to start from a suitable algebra and proceed to find a natural 
generalization of the one species ASEP.We will show that such a natural generalization 
exists. We will obtain the algebra, its representations and many of its properties. In this 
p-species ASEP each particle of type j, hops with rate Vj to its right neighbor site , and 
when it encounters a particle of type i with Vj > Vi they interchanges their site with rate 
Vj — Vi, as if the fast particle overtakes the slow one. In any other case the attempted 
move of the original particle is rejected. The model seems to be relevant as a simplified 
model of one way traffic flow. We consider open boundary conditions where particles ( 
or cars ) enter the system ( the highway ) with a rate proportional to their speed and 
leave at the other side. 

The paper is organized as follows: In section 2 we introduce the algebra and solve 
its'consistency conditions. In section 3 we obtain the Hamiltonian which corresponds 
to this algebra . In section 4 we obtain the uncorrected steady state and finally in 
section 5 we construct the infinite dimensioal representations of the algebra and present 
some of its useful properties. In all the steps we compare our resulsts with the one species 
ASEP. 



2 The Algebra 

We begin with an algebra generated by the elements Di,i = l,...p and E, with the 
following relations: 

D t E = —Di + E (1) 

Vi 

DjDi = ZjiDj + TfciA j > % (2) 

where the parameteres Vi are finite real numbers and the parameters ^ andr/^ are to 
be determined. In writing (1,2) we have had in mind a totally assymetric exclusion 
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process. In order to be consistent with associativity we have to check that 

(DjDi)E = Dj(DiE) (3) 
Using (1) and (2), this requirement determines ^ and as 

& = — Vji = — 4 

Vi - Vj Vi - Vj 

For the present we take all the Vi to be different. Later we will relax this condition. We 
should also check that 

{D k Dj)Di = D k (DjDi) k>j>i (5) 

Using (2) we find that this relation imposes no new relations on the parameters. Since 
any monomials of the form D^Dj...Di with k > j > ...i can be reduced to a linear 
combination of generators, associativity is guaranteed for all the monomials of higher 
degree. The final form of the algebra reads 

DiE = -Di + E (6) 
DjDi = —^—(viDj-VjDi) j>i (7) 

Vi Vj 

where the parameters v\....v p are free. 



3 The Hamiltonian and the Process 

We consider a one dimensional chain of N sites. The Hilbert space of each site is p + 1 
dimensional. The Hamiltonian is written as: 

k=N-l 

H=h 1 + H B + H N = h 1 + £ h* M1 + h N (8) 

k=i 

Where H B is the bulk Hamiltonian and h 1 and h N are boundary terms. In the matrix 
product formalism, the steady state of this Hamiltonian can be written as: 

\P >=< W\V(g)V(g)V(g) ....V®V\V > (9) 
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(E \ 
£>i 



where D 



and \V > and < W\ are vectors in an auxiliary space. According to 



(10) 

(11) 
(12) 



where we take X to be X = 



MPA, we have H\P >— if the following relations are satisfied: 

h B V ®V = X ®V -V® X 

(h N V-X)\V >= 
< W\(h}V + X) = 
/ -1 \ 

with c-numbers. Anticipating the process 

\ x p ) 

from the form of the algebra,we write the bulk Hamiltonian as follows: 

p p 
h B = -J2 Vi( E oi ® Eio ~ E u ® E m ) - J2 yi^j ® E j{ - <g> En) (13) 

i=l j>i 

Here the matrices E^ act on the Hilbert space of one site and have the standard defini- 
tion: (Eij)kj = S^Sji. The first term expresses an ASEP for each species of particles (i) 
with rate t/j. The second term represents an interchange of particles of type (j > i), with 
rate y^-. The parameters yi and i/ij are to be determined. 

The natural choice of the boundary terms should be such that particles of type (i) 
are injected at the left and extracted at the right with their specific rates say a.^ and % 
respectively. So we take: 



h 1 = -J2 a i(E i0 - E 00 ) 



h N = -Y,Ji(Eoi-E u ) 



i=i 



Inserting (13) in (10) leads to the following equations: 



V% 



DjE 



P 



A 
p 



+ XiE 



P 



-DjDi 



x j D 'i x ^ J^) j 



(14) 
(15) 

(16) 
(17) 
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Comparision with (6) and (7) shows that: 



- V A - 

■^i Hi 

V 



(18) 



Inserting the values of the parameters yi and yij in (13) clearly highlights the physical 
process governed by the Hamiltonian. Denoting a vacant site by the symbol d> and a site 
occupied by a particle of type i, by the symbol A iy the process defined by H B is : 



Ad) 



AjAi 



with rate Vi 
AiAj j > i with rate Vj — V; L 

In order for all the rates to be positive we restrict the range of v^s as: 

V 1 <V 2 < V 3 .... < Vp 



(19) 
(20) 



(21) 



Inserting (14) and (15) respectively in (11) and (12) we find the the following explicit 
matrix equations: 



/ -71 -72 -73 

7i 

7 2 . 

73 















"7p 



\ 









/ E \ 



L 

V 



D 



( "I \ 



V 
V 



\V >= 



(22) 



< w\ 



( «i + a 2 + ■■■Oip 

-ai 

-a 2 

-a 3 



V 



-OLr. 



These equations then yield: 







DAV >-- 



\ 






( E \ 

D 1 



( - 1 \ 



+ 



v 

V2 



= 



(23) 



K v f J 



w > 



li 



< W\E =< W\ 



pati 



(24) 
(25) 



6 



and 

* + + = 1 (26) 
V 

Equation (26) which means that the average speed of the particles is unity, is not really a 
condition on the rates, since by a rescaling of time one can always obtain this condition. 
Hereafter we will write this condition simply as < v >— 1. Equation (25) however implies 
that — is independent of i. Thus we set a, = -cm. 
In order that equation (24) be consistent with (7), we require that : 

DjDi - ^ _J - j |^ >= (27) 

After using (24) and rearranging terms, this condition shows that 

li-vi = 7,- - Vj. 

Thus we set 

7i = Vi + P - 1 



Therefore we have 

A|^>= ^ — r|T/> <W\E =<W\- (28) 

+ p — 1 a 

The meaning of the parameters a and j3 are now clear. From (26) we find: 

7i +72 + -7 P fl 

P 

Thus a and /3 are respectively the total rate of injection and the average rate of extraction 
of particles. In order for the individual rates to be positive we restrict further the range 
of speeds as 

1 -P < v, < v 2 < v 3 .... < v p (29) 

The above results clearly justifies this model as a natural generalization of the 1-species 
ASEP. 

We conclude this section with formulas for the current operators. A simple way to 
obtain these, is to directly refer to the process (19-20) and write the equation for the 
one point function < rj^ > which is the average density of particles of type (i) at site k. 
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Determining from (19-20) various ways of increasing and decreasing of this density and 
denoting the probabilty of site k to be vacant by we find: 

j t < r k ] > = Vi< r^efc > -Vi < 4 l) e k+1 > 

j<i j<i 

~ Efe - «0 < t&t® > + E(^ - Vi) < ri j) rtl > (30) 

j>i j>i 

Note that inside any correclation function we have = 1 — J2i T k ' > - Equation (30) can 
be rewritten as a continuity equation 

± < r (0 >_< j(<) > _ < j(<) > 

where 

< Jf >= «i < r&e* > + E(^ - «i) < > " E(«i " «*) < > ( 31 ) 

j<i j>i 



According to MPA 



< 4" >= 



< w\c k - 2 j®c N ~ k \v > 

< W\C N \V > 



where C := E + := E + ^j(A + A + ■■■D p ). The current operators can be read from 
(31) to be: 

J W = v^E + J> - ^)^ - E(«i - — ~ (32) 

p p p U vp 

Using eqs.(6,7), we obtain: 

jw = ^(1 A + E)~h E(^A - ^A) - ^ EKA - «<A-) 
p y » p j<i f j>j 

= -a + -^-^(E^)a + ^(E^) 
p p p 2 f^i p u 

= -Di + -ViE - \(p - Vi)Di + ^(^ - A) 
p p p z p z p A 

= -C (33) 
P 
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We thus obtain 



< J (i) >-- 



v i <W\C N - 1 \V > 



(34) 



p <W\C N \V > 

Therefore all the currents are simply proportional to the average current J, however J 
itself has a highly nontrivial dependence on the hopping rates. 

In the next section we will find the one dimensional representations of the algebra. This 
case corresponds to the steady state being characterized by a Bernouli measure. Although 
very simple this steady state has a rather rich structure. 



4 One Dimensional Representations and the Uncor- 
related Steady State 

In the one dimensional representation the operators Di and E are represented by c- 
numbers A, and e respectively. From equation (28) we have : 



A; 



Vi + P-1 



1 

a 



(35) 



combining (35) and (6) gives the condition a + (3 = 1 on the average rates which is of 
the same form as in the one species asep. This, together with (29) means that for an 
uncorrelated steady state to exist, the minimum speed v\ should be greater than the 
average arrival rate a. The steady state is now given by \P >= \p >® N where 

/ e \ 

Al 

L 

V 



\p>-- 



c = e+ -(Ai + A 2 ...A P ) = e + -A 
p ' p 



(36) 



V f J 

The density and current of particles of type (i) and the total density and the total current 
are all site independent and are respectively given by (see 34) 



p(a,i) 



V 

v 



J(a, i) 



p 



p(ot) 



A 
P 



p 



J{a) 



p 



(37) 
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It is better to consider the limiting case p — > oo, that is, we assume that particle 
speeds are taken from a continous probability distribution P(v). Condition (26) is then 
transformed into 

< v >: = / vP(v)dv = 1 

Discrete quantities ^f(i) are transformed into the continous functions f(v)P(v) and sums 
into integrals. Instead of (36) we have 

A(a,v)P(v) T . N vP(v) , \ 

p(a, v) = v ' \ ! J(a, v) = 38 

IK 1 e + A(a) v 1 e + A(a) v 1 

p(a) = A( A ^ J(a) = 1 —— (39) 

IK 1 e + A(a) v 1 e + A(a) v 1 

where 

A(a,v) = ^^ and A(a) = f^^P(v)dv (40) 

v — a J v — a 

As an example, in the following,we take a distribution which at low speeds vanishes as 

some power of v — a, so that a is indeed the minimum speed of particles. Requiring that 

the distribution has an exponential decay rate A at high speeds and the average speed 

be unity gives: 

P\(v) = : t; 7 - a^e^ 1 (41) 

; A m+1 r(m + l) v ; v ; 

where m + 1 = ^^.Here m is not necessarily an integer. A similar choice of P(v) has 
also been made in [14]. Note that since m is to be positive for each choice of A, we should 
have < a < 1 — A. This distribution is peaked at v — 1 — A. Inserting (41) into (40) 
and doing the integral gives: 

A(a) = - - (42) 

A + a — 1 

from which we obtain 

a(A + a-l) a(A-l) ,. oN 

A — 1 + aA A — 1 + aA 

Eliminating a between the above equations gives the current versus density: 

J w = "( 1 -T^rt) (44) 

The curves J(a) and J{p) are shown in figs.(l and 2) for various values of the widths of 
the distributions. Note that for zero width ( A = 0) the familiar results of the one species 
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ASEP are obtained, namely: p = a and J = p(l — p). For any finite value of A p is an 
increasing function of a but J (a) has a maximum at a := ^y^(l — y/l — A) The current 
at this optimal value of the arrival rate is 

^(A) = (l-A)^ 1 "^ 1 ) (45) 

The maximum of J as a function of p occurs at p — — A — (1 — A)). It is seen 

that the maximum current has it is highest value of \ only when all the particles have 
the same hopping rates. It is interesting to note that even in such a simple uncorrleated 
steady state, a variance in hopping rates reduces the current. 

The final quantity we consider is the average density of particles of different speeds, p(a, v) 
as a function of v. From (38-41) we have 

p(a, v) oc v(v - a) m ~ l e ^ ~ (46) 

The interesting point is that when m > 1 ( or a < 1 - 2A) the density vanishes at the 
lowest speed a and has a maximum at v max := |(1 — A + \J(1 — A) 2 — 4Aa). However 
when m < 1 or (a > 1 — 2A) the distribution changes and the density of the slowest 
particles rises abroptly . This transition is one of the interesting features of this process 
and is somehow reminicent of the Bose condensation first noted in [14]. 



5 Representations 

As in the one-species ASEP the representations of this algebra are either one dimensional 
or infinite dimensional. To see this we first note that there is no non-zero eigenstate 
of £",say |e > such taht E\e >= > for any %. Since if there is any such vector 
then by acting on both sides with D-i and using (5) we obtain —\e >= and since 
Vi is finite this means that |e >= 0. Following [2] we know show that in any finite 
dimensional representation of the algebra all the generators mutually commute. For if 
the representation is finite dimensional then the matrices E — ...E — ^- having no 

zero eigenvalue, are all invertible and hence (5) gives Di = E(E — — ) _1 which in turn 
means that the representation is commutative. Thus the representations of this algebra 
are either one dimensional or infinite dimensional. Note that in the one dimensioanl case 
the expression Di = E(E — ^-) _1 automatically satisfy (6,7). Thus such one dimensional 
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representations really exist, and were used in section (4) to find the uncorrected steady 
state. 

To find the infinite dimensional representations we assume that there is one single vector 
denoted by |0 > such that 



DAO >= dAO > 



1, ...p 



(47) 



where the paramters di are to be determined. We then consider the vector space W 
spanned by the formal vectors {\n >:= E n \0 >,n = 0, 1,2, ....}. Clearly 

E\n >= \n+ 1 >, Vn (48) 

Iterating (5) we find 

DiE n = vr n Di + v~ n+1 E + v^ n+2 E 2 + ...v^E'^ 1 + E n (49) 

Thus we obtain: 

Di\n >= v- n di\0 > +v^ n+1 \l > +vr n+2 \2 > +...v^\n - 1 > +\n > (50) 

In order to chech eq.(6) in this representation it is enough to consider only the state |0 >, 
since all the other states are obtained by acting on this state using the algebraic relations 
which have been found to be consistent. Therefore we require that 



DjD, t 



Vi - v, 



{viDj -VjDiWO >= 



(51) 



this fixes the parameters di to be di = where e is a constant independent of i . The 
explicit matrix form of the generators are then as follows: 





(0 . . . . 


\ 




10... 






.10.. 




E = 


..10. 






v 


• • / 








Di 



di. 
1 



di 


di 


di 






~t 




1 


1 








1 


{_ 

Vi 


i 



(52) 



Remarks: 
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1) If e = 0,then all the d{ — 1 and using (50) one sees that DiDj = DjDi Vi, j al- 
though DiE 7^ EDi. Hereafter we restrict ourselves this case. Note that this does not 
mean that the two point functions < r^V/^ > and < t^t® > are equal. 

2) When e = we have Dj = Dj for = Vj, so that equation (7) becomes vaccu- 
ous in this case and no singularity arises due to Vi = Vj in this equation. Furthermore 
one can safely eliminate one of these generators for the other one. This means that the 
(p-l)-ASEP algebra is naturally embedded in the (p)-ASEP algebra. 

The ket and bra vectors \V > and < W\ are found to be: 

oo 

|V>=£(1 -/?)>> (53) 

n=0 

oo 

< w\ = J2( a )~ n < n \ ( 54 ) 

n=0 

Having a representation at hand is only half of the way in obtaining explicit expres- 
sions for the physical quantites like currents and correlation functions. To obtain these 
quantities,one should do lenghty calculations on the matrix elements of porduct of op- 
erators. Our work is not complete in the sence that we do not yet have the final explicit 
expression for these matrix elements. However we explore the properties of the represen- 
tations as far as we can, hoping that using these, the rest of the problem be solved in 
another occasion, by this author or by others. 

5.1 Eigenvectors of Di 's 

Let z be any complex number, with \z\ < V\. Define: 

oo 

\z >= z n \ n > (55) 

?1=0 

Then \z > is a common eigenvector of all ZVs: 

Di\z>=-^—\z> (56) 
Vi - z 

This can be easily proved by direct calculation. The states \z > can in fact be thought of 
as the coherent states of the algebra constructed as follows: \z >:= |0 > where 
is understood by its power series. We now define the dual vector < z\ — J2^=o z ~ n < n l- 
Then we have the following theorem, the proof of which is accomplished by noting that 
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< n\E =< n — 1| and < 0\E — 0, and doing straightforward calculations. 
Theorem: 

a) < z\E — - < z\ 

b) The coherent states form a basis of the space W, and: 

\n >= [ —dzdzz~ n \z > 

J\z\<l 7T 

c) The coherent states \z > form an over complete basis for W. A complete basis is 
obtained by taking only the states with fixed \z\, with a completeness relation: 

dz . 

\z X z\ 



fixed\z\ 2niz 



d) < uo\z >= jzt for \z\ < \oj\ 



Remarks: 

1) It is now appropriate to recall that the states \V > and < W\ are in fact coherent 
states and so it is best to denote them respectively by |1 — (3 > and< a\. 

2) The calculation of any physical quantity is now reduced to the calculation of matrix 
elements of powers of C between the coherent states. For example, the average density 
< n\ > is given as: 

<n i >= J_/ Jt_ <a \c k - 1 \ z><z \C N - k \l-(3>-^- (57) 



Zn J(l-P<\z\<a, \z\<Vi) 2niz ' ' Vi — Z 

where Zn =< a\C N \l — (3 >. Or the probability of a segment [k,l] be empty is given by 
P(fc,Z):=^-/ < «|C fc_1 k >< z\C N - k ~ l \l-P > -^r (58) 



Z N J\-p<\z\<a, 2niz 1 z l k 

Once the eigenvalues and eigenvectors of C are known, the problem can be solved 
completely. In the following we will give an implicit formula for these objects. 
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5.2 Eigenvalues of C 

First we consider the case p — 1. For any coherent state \z > we have E\z >= -(\z > 
— 10 >). For any two complex numbers z and uj define: 

\z,uj >= z\z > —uj\uj > (59) 

Then it follows that: 

E\z,uj >= \z > -\oj > (60) 
Since for p — 1 we have D\z >= -r—\ z > we obtain: 

1 1 



c 



z,u >= (D + E)\z,u >= \z>-- la; > (61) 

l — z 1 — UJ 



Therefore the state \z, oj > is an eigenstate of C if we have z(l — z) — ui(l — uj), the 
solution of which is u = z or u; = 1 — z The second solution is acceptable and hence we 
have: 

C\z,l- z >= — -\z,l-z> (62) 

z{l — z) 

Now let p > 1. Using (56) and (60), we obtain: 

C\z,u >= zt)(z)\z > — wr](uj)\uj > (63) 

where 

1 1 v v- 

^) = - + -E— (64) 

The state \z, ui > is an eigenvector of C if z and w lie on the curve r](z) = r)(u>). We will 
then have: C\z,u(z) >= r)(z)\z,u>(z) >. 

Note that in the p — ► oo z and uj are related by the following equation: 

i + fvm = i + rvm (65) 

z J v — Z UJ J v — UJ 

In order to understand the connection with the eigenstate obtained in Derrida et al 
[2] in the p — 1 case we change the basis of W as: 

\e n >= (E — l) n |0 >, | e „>=|0> (66) 

Using ( 6,7) it is easily found that: 

E\e n >= \e n > +\e n+1 > 
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D\e n >= \e n > +|e n _i > 
C\e n >= |e n _i > +2|e„ > +\e n+1 > 



Thus we have 



1 . 1 . 1 / 1 \ , 1 . z 

\z >= >= ; dO >= >= > d 

1 1-zE 1 I — z — z(E — 1) 1- z\l- j^EJ 1 l-z l l-z 

(67) 

where by \z >d we mean \z >a— Y^=o zn \ e n > Therefore 



1 — 2 1 — z z z n=0 V 1 — z z J 

(68) 

Taking = e l6 gives 

oo 

\z,l- z >=^2sin(n + l)0\e n > (69) 

n=0 

which is denoted by \6 > in [2] with eigenvalue z ^_ z ^ = 2(1 + cosQ). This concludes our 
treatment of the representations of the algebra. 

6 Discussion 

This work can be pursued further in the following directions: 

a) Finding a solution of the equation rj(z) = t](uj) either for a low value of p {e.g.p = 2) or 
in the p — >■ oo limit, and then expanding the state |1 — f3 > in terms of the eigenstates 
of C. In this way one will obtain the current and the phase diagram of the system. 

b) Finding a solution of the mean field equations, either numerically or analytically. 
These equations can be written down from (31). In the p — ► oo limit they are: 

avei = J k (v) = J k +i(v) = ... = (v + (3 - l)n N (v) 

where 

J k (v) = vn k {y)e k+1 + / (v - v')n k (v)n k+1 (v')dv' - / (v' - v)n k (v')n k+ i(v)dv' 

JO Jv 

Here n k (v) is average density of particles of speeds v at site k and e k := 1 — / °° n k (v)dv 
is the probability of site k being vacant. In the one species model it is known that the 
mean field analysis gives the phase diagram correctly. So it will be interesting to see for 
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a typical probability distribution how the phase diagram will be modified. 
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Figure Captions 



Figure 1: The current versus the arrival rate of particles for different values of A. 
Figure 2: The current versus the density for different widths of the distribution. 
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